Scaling of the irreducible SO(3)-invariants of velocity correlations in turbulence 
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The scaling behavior of the S0(3) irreducible amplitudes 
d„{r) of velocity structure tensors (see L'vov, Podivilov, and 
Procaccia, Phys. Rev. Lett. (1997)) is numerically examined 
for Navier-Stokes turbulence. Here, / characterizes the irre- 
ducible representation by the index of the corresponding Leg- 
endre polynomial, and n denotes the tensorial rank, i.e., the 
order of the moment. For moments of different order n but 
with the same representation index I extended self similar- 
ity (ESS) towards large scales is found. Intermittency seems 
to increase with I. We estimate that a crossover behavior 
between different inertial subrange scaling regimes in the lon- 
gitudinal and transversal structure functions will hardly be 
detectable for achievable Reynolds numbers. 



The most fundamental objects to analyze the struc- 
ture of turbulent velocity fields u{x,t) are the tenso- 
rial moments of the velocity differences Vi{r;x,t) = 
Ui{x+r, t) — Ui{x, t), averaged over time t or/and position 
X, considered as functions of scale r, 



(1) 



If the eddy size r = |r| is in the inertial subrange (ISR), 
i.e., T] r <^ L, algebraical scaling of the moments is 
expected. Here, rj is the inner (Kolmogorov) scale and L 
the external length scale 0|. If the turbulent flow field 
can be considered as statistically isotropic (or close to), 
one better uses rotational invariants instead of the ten- 
sorial components, in order to cope with the multitude 
of scaling exponents. The most commonly used invari- 
ants are the structure functions of the longitudinal veloc- 
ity component vl = f • r° and the transversal velocity 
Vt = V — ult"; here, denotes the unit vector in r 
direction. We denote these structure functions as 



advances in experimental technology (^-|^ and computa- 
tional power and technique |^ raised increasing doubts 
if this is true for general moments of order n, as it is for 
the most often considered 2nd order structure function, 
n — 2, where the condition of incompressibility enforces 
D2 oc D2 oc (X r'»^. For general n, it was found 
in several experiments and simulations that the degree 
of intermittency (i.e., the deviations of the scaling expo- 
nents from the classical value Cn = 'T'/S) is considerably 
larger in the transversal moments compared to the lon- 
gitudinal ones; for a summary of the results see table 1 
ofref. §. 

In a recent paper, L'vov, Podivilov, and Procaccia 
Q suggested that it were not the longitudinal and the 
transversal structure functions which obey clean alge- 
braic scaling, but the amplitudes of the moment tensor 
eq. (0) decomposed into the irreducible representations 
of the rotation group SO (3), 



d^r) oc ((t;2(r;a;,t))"/2p;(^o . ^O)^ ^ 
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The representation label I runs through < I < n with 
the same parity as n, if statistical reflection symmetry 
of the turbulent flow fleld is guaranteed ||]; Pi is the 
Legendre polynomial. The amplitude of the unity repre- 
sentation, d°(r), is already part of the conventional set 
of structure function, since d^{r) oc {r). 

For the second and fourth order structure tensors the 
amplitudes d\{r) and d\{r) are linear combinations of the 
longitudinal, transversal, and modulus structure func- 
tions. We follow L'vov et al.'s definitions Q ao = dj, 
ci2 — d^, C{) ~ d^, C2 — d\, Cj^ ^ d\ obtaining 
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Di:{r)^{[vL{r;x,t)r)^r'^^ 



D^{r) = {\vT{r;x,t)n^r'^^; 



(2) 



(3) 



both are assumed to scale in the ISR with the correspond- 
ing exponents and (J. A third convenient structure 
function is the n-th order moment of the modulus of the 
eddy velocity difference v(r;x,t) which again is assumed 
to scale 



D^'{r) = {\v{r;x,t)n^r<"^ 



(4) 



Traditionally, it was believed that all three scaling expo- 
nents are the same, (n = Cn = Cn = • But recent 
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On the rhs also other ways of representing the n-rank 
velocity correlation tensor can alternately be given, us- 
ing e.g. D\ in (^ or mixed transversal/longitudinal mo- 
ments in (|^) as done in eq. (13.81) of ref. ||l|] or in reL H 



which uses Dn 
-D1122 = {vfvj) 



D^, D22 



2 /2, i^iiii = {Vli - , 



and D2222 = 



= 3D 
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where the 1-axis has been put in the longitudinal direc- 
tion parallel to r. For these structure functions we obtain 
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The point of L'vov et al. is that the invariants ai, ci on 
the Ihs of eqs. (||), (Q), and (^) are distinguished because 
the djj(r) are the ampHtudcs of the structure tensor for 
its decomposition into the components of the irreducible 
representations of the rotational symmetry group SO (3). 

In this paper we present the scahng properties of the 
fourth order moments d^^ir) from a fuU numerical simu- 
lation of the Navier-Stokes equation on a 96^ grid with 
periodic boundary conditions. The numerical turbulence 
is forced on the largest scales, the averaging time is about 
120 large eddy turnovers, and the Taylor-Reynolds num- 
ber is Re\ = 110. The isotropy of the flow has carefully 
been checked; for details of the simulation we refer to ref. 
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FIG. 1. Fourth order structure functions co,2,4(»') as func- 
tions of r. As tested for smaller Reynolds numbers and longer 
averaging times, the wiggle in C4,{r) at large r is not statisti- 
cally safe. It seems that very long averaging times are neces- 
sary for moments with large I to converge at large scales. The 
inset shows co,2,4(r') (top to bottom) as they follow from the 
Batchelor parametrization and eqs. (2) of ref. . The dif- 
ferent magnitudes and the distinct transitional behavior of the 
different irreducible representations can be recognized. Note 
that the local slope of lgC4(r) vs Igr around the transition is 
not monotonous. - By assumption the ISR scaling exponents 
are the same. 

The second order moments all asymptotically scale 
the same because of incompressibility. Assuming clas- 
sical scaling (2 = 2/3 one obtains D2 — 4L)J/3 and 
Go — lla2 = I?2^/3. In figure |l| we give the fourth or- 
der structure functions co,2.4('")- As expected for this low 
Rex, the scaling properties of these structure functions is 
very poor, because there is not yet a well developed ISR. 
There is analytical behavior oc r"* in the viscous sub- 
range (VSR) followed by a transition and leveling off in 
the inertial and stirring subrange around r ^ L. What 



can be said, however, is that with increasing / (i) the 
magnitude of c; (r) decreases and (ii) the degree of inter- 
mittency seems to increase, < Cl < C4 < 4/3. The 
reason for (i) is that cq is a sum of positive definite struc- 
ture functions, whereas C2 and the more C4 are differ- 
ences thereof, similar to 02 = (-D22 — ^ii)/3 which also 
is much smaller than oq = Dn -\- 2D22- The reason for 
(ii) presumably is that larger / in (^) means the probing 
of smaller scale structures which are traditionally associ- 
ated with stronger intermittency. 

Fortunately, the extended self similarity method (ESS, 
|l0| ) allows for more quantitative statements. Here, we 
focus on the scaling of the fourth order structure func- 
tions vs second order ones. More specifically, to visualize 
the deviations from classical scaling we calculate com- 
pensated ESS plots D\/{D\f vs D\ 



^4/(^2)^ vs d\, 1 — 0,2, see figure I 



i — L,T, M, and 
For I — we find 
ESS scaling from r ~ lOry up to r ~ L, resembling the 
ESS scaling for the longitudinal and transversal struc- 
ture functions figure ||a which was extensively analyzed 
in |l^, For / = 2 we find ESS towards large scales 
r > 5O77, but no ESS towards smaller scales r < 50r]. 
Instead, there is a bump in the curve 02/03 vs 02 for 
r ~ 3677. We checked very carefully whether the bump 
would smooth for increasing averaging time. This is not 
the case. It also persists for different type of large scale 
forcing and smaller Reynolds number, but much larger 
averaging time. 

At first sight, the bump was a surprise to us. How- 
ever, we suggest that it can be understood as a transi- 
tion phenomenon from the VSR to the ISR, similar to 
the one seen in figure 3 of ref. |ll]. Hitherto, it was 
not observed in ESS plots of longitudinal and transversal 
structure functions as both are dominanted by the rank 
zero contribution which does show ESS. 

For further support of this interpretation, we 
parametrize the d'^{r) within Batchelor's parametrization 

i, 



d"„{r) (X D^'ir) (X r" 



(-«+C„)/2 
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n = 2, 4, with the She-Leveque model ||l2| values (2 = 
0.70, C4 = 1-28, and rc = lOry. Incompressibility gives 
Z52 (r) = r~^ D^{r)Pdf and via eq. ^) all other sec- 
ond order structure functions follow. For the fourth order 
moments, an analogous relation does not exist. However, 
within some closure approximations |^ (whose nature 
is controversal) , all 4th order structure functions follow 
from D4 (r) , cf . eq. (2) of ref. |^ . We stress that those re- 
lations are not generally true and their consequence that 
all 4th order structure functions scale the same is in di- 
rect contradiction to our and others' findings. However, 
for the demonstration of transitional effects, for which 
the different intermittency in the ISR does not matter, 
eqs. (2a) - (2c) of ref. Q could be useful. Employing 
them we derive an ODE for Dim = -Df, 
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FIG. 2. Compensated ESS plots for different space direc- 
tions of fourth vs second order structure functions, giving the 
ISR scaUng exponent pi = {Q — 2C,i)/C2- The plot C2/a\ vs 02 
displays a distinct bump around 35r; for all three space direc- 
tions, though it is slightly differently developed in strength, 
possibly because of the unavoidable anisotropy in the forcing, 
possibly because of the too short averaging time (Fig. 2b, up- 
per). - Convergence is much less a problem in the ESS plot 
co/flo '^s ao (Fig. 2b, lower) and consequently also in the ESS 
plots of the longitudinal and transversal structure functions 
where the rank zero contributions dominate (Fig. 2a). 



(10) 



which can numerically be solved for the given D^^ (r) of 
eq. (pD. With 



D1122 = -D^ + ——D^, 



(11) 



Do 



Di 



—r—Dl + —r^—D^ (12) 
16 ar 16 ar^ 

and eqs. (^ - (||) all other fourth order moments follow. 
The resulting compensated ESS plots are shown in figure 
^. Indeed, in the ESS plot C2/a\ vs 02 a bump occurs 
at the VSR-ISR transition, similar to what we found in 
the numerical simulation. All other shown compensated 
ESS plots are dominated by D^' — \/5co and I?!^ = 3ao, 
and therefore display ESS, as cq vs ag does. After the 



functional shape of the ESS transition from VSR to ISR 
in Z?4^ vs is now believed to be rather universal, it 
would be worth while to analyze in various experimental 
and numerical flows, whether the first angular contribu- 
tion, i.e., C2 vs 02, also is somehow universal and thus 
displays the type of structure as we found in figure ^b. 

We now come back to the numerical results and focus 
on the ESS scaling exponents of figure || which we denote 
by P^ = (CI - Kl)/C,l I = L,T,M or i = 0,2. The 
deviation of the pi from zero characterizes the degree 
of intermittency of the corresponding moment. We find 
PL = -0.15, pt = -0.30 and pM = Po = -0.25, p2 = 
—0.5, again showing that the degree of intermittency is 
higher in the d}^ with larger I. The She-Leveque model 
value (with the original She-Leveque parameters adopted 
to the longitudinal structure function) ||l^ for p is p = 
-0.16. 

We checked the possibility of scaling behavior if ampli- 
tudes corresponding to different irreducible subspaces are 
mixed: We do not find ESS if we plot structure functions 
1^4/(^2)^ vs c?2 with different I ^ I'. 

It will not have escaped the reader's attention that the 
simultaneous assumption of pure scaling behavior of both 

L T AI 2 4 

the Da ' as well as the c^' ' is self contradictory if the 
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^4 as well as the C4 

exponents with different I are different. We follow L'vov 



et al.'s argument that the d„ are the more fundamental 
structure functions and (employing eq. (Q) and incom- 
pressibility) write the ratios ' 
ratios of the di, , 



/{D^'^^Y a sum of 
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From our numerics (see fig. ^ the first term is found to be 
the leading order term. It represents the scaling of the 
modulus structure functions eq. (^. In the first (and 
larger) correction term the approximation ao ~ lla2 
(resulting from C2 ~ 2/3 and incompressibility) can be 
made, leaving only ratios whose scaling we can deter- 
mine from the ESS-plots figure |^b; (the C4-term hardly 
contributes for large r). With that approximation the 
qualitative features of figure ^ja can be understood from 
eqs. ( |l3|) and (|l4|): As C2(r) < the 02/03 correction 
term to the leading cq/gq term is negative [positive] for 
Di/iD^)'^ [DJ /{D'^)'^], leading to a less steep [steeper] 
"apparent" slope for the ESS exponents pL — —0.15 
[pt = —0.30] of the longitudinal [transversal] structure 
function compared to the leading contribution with pQ — 
Pm = —0.25. Even that the correction to po = —0.25 is 
twice as big for the longitudinal structure function than 
for the transversal one can be seen from eqs. (|l3|) and 

Finally, we would like to estimate for what Reynolds 
number two distinct scaling regimes (in r) may be ob- 
servable in ' . Therefore, we plug in the scaling laws 
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FIG. 3. Compensated ESS plots following from the Batch- 
elor parametrizations (|^) and eqs. (2) of ref. j^. The curve 
02/02 vs a2 is the only ESS plot which is not dominanted 
by the tensors of rank zero and it does show a bump at the 
VSR-ISR transition. That the ISR scaling exponent pi is the 
same for the four curves shown here for demonstrative reasons 
is a trivial consequence of the closure assumptions of ref. 
but does not hold for the real data. 



C2Pi 



0,2 
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and obtain with the numerical values at r = L, cq/oq « 6 



and C2/a\ 



'100, 



r \ C2(p2-pq) 



C2PT 4 

fX FFT-TT oc r 



1 / 7- \ C2(p2-po) 



C4-corr.J , 

(16) 
+ C4-corr 

(17) 



We get a - 0.2, (2(^2 - Po) = §(-0.5 + 0.25) w -0.17. 
Note that for small enough r the second term in ( |l7| ) 
may dominate the first one and for even smaller r the 
third term will contribute. [In eq. (|l^) the situation is 
more complicated as the second term has negative sign, 
but the Ihs is positive definite.] Therefore, in princi- 
ple D\I(I)\Y shows several different scaling regimes. 
However, it will be very hard to detect these different 
regimes as the required span of the Reynolds number is 
too large. In eq. ( p7| ) the ratio Ljr has to be as large as 
Ljr = (2/q;)^/°-^^ ^ 10^ for the second term to overtake 
the first one. We put r = -q and estimate that this means 



Re ^ 10*, which is hard to achieve in today's experimen- 
tal or numerical flows. What shall be detectable if L'vov 
et al.'s conjecture is right is that the apparent scahng 
exponents of the structure functions D^'^{r) or ESS scal- 
ing exponents thereof are slightly Re dependent whereas 
the scaling exponents of the irreducible objects (i^(r) or 
their ESS exponents (the exponents of plots |djj(r)| vs 
|(ij,j(r)|) might well be universal, i.e., Reynolds number 
independent. Measuring such exponents up to very large 
Reynolds numbers j^] seems to be of prime importance 
for our further understanding of fully developed turbu- 
lence. 
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